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Abstract 

Linearized equations of motion for gravitational and scalar fields are found and solved 
in a stabilized brane world model in five-dimensional Brans-Dicke theory. The physical 
degrees of freedom are isolated, the mass spectrum of Kaluza-Klein excitations is found 
and the coupling constants of these excitations to matter on the negative tension brane 
are calculated. 

1 Introduction 

Nowadays models with extra dimensions with the fundamental energy scale lying in the TeV 
range are widely discussed in scientific literature. One of the most known and interesting from 
the phenomenological point of view is the Randall- Sundrum model pQ. It describes two branes 
interacting with gravity in a five-dimensional space-time and provides an original solution to 
the hierarchy problem of gravitational interaction [H[2j[3]. Nevertheless, the Randall- Sundrum 
model possesses an essential flaw: the distance between the branes is not fixed by the parameters 
of the model. This leads to the existence of the massless scalar field - the radion - in the four- 
dimensional effective theory on the branes. The coupling constant of this field to matter on 
the negative tension brane, which is assumed to trap the Standard Model fields, appears to be 
very large, which contradicts experimental data even at the level of classical experiments [21 II]- 
This problem was solved by introducing an extra scalar field living in the bulk. The most 
consistent model was proposed in paper [5], where exact solutions to equations of motion for 
the background metric and scalar field were found. The size of the extra dimension is defined 
by the boundary conditions on the branes. Nevertheless, the background solution for the metric 
in this model is more complicated compared to the simple solution in the Randall-Sundrum 
model. It turned out that the corresponding linearized equations of motion in this model 
can be solved analytically only for a certain choice of the model parameters, i.e. when the 
background solution for the metric can be approximated by the solution of the unstabilized 
Randall-Sundrum model [BJ. A question arises: is there a stabilized solution in the system of 
two branes admitting the simple Randall-Sundrum solution for the metric? In paper [7] it was 
found that one can obtain such a solution in the case of non-minimal coupling of stabilizing 
scalar field to gravity. But the background solution for the scalar field in this paper has a rather 
complicated form. 

One of the standard forms of the non-minimal coupling of scalar field to gravity is the linear 
interaction with the scalar curvature used in Brans-Dicke theory. It has been found that it 
is possible to get a stabilized model with two branes in five-dimensional Brans-Dicke theory, 
admitting the simple Randall-Sundrum solution for the metric, whereas the solution for the 
scalar field also has the form of a simple exponential function [8]. 



In the present paper we study linearized gravity in the stabilized brane world model proposed 
in [8]. It turned out that due to the simplicity of the background solution, the linearized 
equations of motion can be solved analytically for any physically interesting choice of the model 
parameters, contrary to the case of stabilized Randall-Sundrum model [5] . We also calculate the 
coupling constants of physical degrees of freedom to matter on the negative tension brane, where 
the hierarchy problem of gravitational interaction is solved, and describe the mass spectrum of 
Kaluza-Klein modes. 



2 Background solution 

Let us consider gravity in a five- dimensional space-time E = M4 x S 1 /Z2, interacting with 
two branes and with the scalar field <fi. Let us denote coordinates in E by {x M } = 
M = 0, 1, 2, 3, 4, where {x^}, fi = 0, 1, 2, 3 are four- dimensional coordinates and the coordinate 
y = x 4 , —L < y < L, corresponds to the extra dimension. The extra dimension forms 
the orbifold S /Z2, which is a circle of diameter 2L/tt with the points y and — y identified. 
Correspondingly, the metric gMN arid the scalar field satisfy the orbifold symmetry conditions 

g^(x,-y) =g tlv (x,y), g^(x,-y) = -g^(x,y), (1) 
g M (x, -y) = gu(x, y), -y) = y). 

The branes are located at the fixed points of the orbifold y = y = L. 
The action of the model has the form 

•L 



S = I d x j dy^f^g 



u g MN d M <pd N <p - v{<p) 



(2) 



g\i((J))d 4: x — / a/— g\2{4>)d A x. 
y=o Jy=L 



Here V(<fi) is the scalar field potential in five-dimensional space-time, Ai i 2(0) are scalar 
field potentials on the branes, u is the five-dimensional Brans-Dicke parameter, g^ v denotes an 
induced metric on the branes. The signature of the metric guN is chosen to be (— , +, +, +, +). 
Subscripts 1 and 2 label the branes. 

We consider the following standard form of the background metric 

ds 2 = lMN dx M dx N = e 2aiy) Vflu dx^dx u + dy 2 (3) 

with r]^ being the flat Minkowski metric, which preserves the Poincare invariance in four- 
dimensional flat space-time, and the following form of the background solution for the scalar 
field 

<j)(x,y) = <j)(y). (4) 
Functions cr(y),(f)(y) are defined by the equations of motion [8j 



" ' " + 4a>') - 4a" - lO^') 2 - (5) 



u ,,/s2 ^ dV 1 d\i ld\2 



<%)-n-7X%-£) = 0- 



20 2Vr ' 2d4> 2 d<p 2 

6(a') 2 0-^^(0') 2 -^)+4aV = 0, (6) 

3a" ( f> + ^(cf>') 2 + ^-a'<P' + ^X 1 S(y) + ^X25(y-L) = 0. (7) 
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Here and below ' = <9 4 = d/dy. 

Let us consider the scalar field potentials to be 

V{cj>) = A<p, Ai, 2 = ±\<f>. (8) 

In this case functions cr(y), <p(y), which are solution to equations @-([7j), take the form 

a = -k\y\+C, <p = C 1 e~ ulyl , (9) 

where 



it A 



(3w+4)(4uj+5) ' 

k = (w + l)u, (10) 

Background solution for the metric has the same form as that in the unstabilized Randall- 
Sundrum model [1]. We consider the case when the matter is located on the brane at y — L, 
like in the Randall-Sundrum model. To this end we take C = kL, which makes the four- 
dimensional coordinates on this brane Galilean (see [21 El II])- 

To fix the size of the extra dimension, let us add the following terms to the brane potentials: 

AA li2 = ^(0-^ 2 ) 2 . (11) 
The equations of motion appear to be satisfied provided 

<fr\ y =o = Ci= Vx, (f)\ y= L = v 2 - (12) 

Thus, the distance between the branes is defined by the boundary conditions for the field 
and can be expressed through the parameters of the potentials as 

u \v 2 J 

We suppose that V\ ~ v 2 and uL < 1. We also suppose that kL pa 35, which can be achieved 
if u > 35. Note that the mechanism of stabilization is based on the dependence of the scalar 
field background solution on the coordinate of the extra dimension. Parameters v\ )2 of the 
potentials, made dimensionless by the fundamental five- dimensional energy scale A, should be 
positive values of the order of 0(1), i.e. there should not be any hierarchical difference. 

As it was noted above, the five- dimensional scalar field usually minimally couples to gravity 
in stabilized brane world models [5]. In the case of non- minimal coupling, there exists a con- 
formal transformation transforming the action from the Jordan frame to the Einstein frame, in 
which scalar field minimally couples to gravity. For the action of form (j2J) these transformations 
are presented in [S] (with the coordinate transformations bringing the background metric to 
the standard form). The difference between the initial and the transformed actions is charac- 
terized by the interaction of matter on the branes with gravity, i.e. by the metric which we 
are supposed to perceive. It appears that if we live in the world, in which five-dimensional 
scalar field is non-minimally coupled to five- dimensional curvature, it is more convenient to 
examine linearized gravity using the initial, untransformed action. It allows one to simplify the 
derivation of the Kaluza-Klein mode mass spectrum and the coupling constants of these modes 
to matter on the branes. But the use of the conformal transformations allows one to simplify 
the choice of gauge conditions, which are necessary for isolating the physical degrees of freedom 
of the theory. We will discuss this point in the next section. 
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3 Linearized equations of motion and the choice of gauge 
conditions 



To study linearized gravity one should derive the second variation Lagrangian of the model. To 
this end let us parameterize the metric and the scalar field as 



9mn(x, y) = 1mn{v) + h MN {x, y), <j>(x, y) = <j>{y) + f(x, y), 



(14) 



where (f)(y) is the background solution of the scalar field, substitute it into action and 
retain the terms of the second order in fluctuations (below we will use the short notation 
for the background solution 4>(y)). The corresponding second variation Lagrangian appears to 
be extremely large and we do not present it here in the explicit form. We only present the 
linearized equations of motion for fluctuations of metric and stabilizing scalar field, which follow 
from this Lagrangian: 



1. /^-component 



[dad^h^ - dyfPhm, - d u d a h afl + d^h^) + (f)d ll d u h+ 
+<pd,j,d v hu - 4>d 4 {dfji 4v + dyhty) - 2a'<f)(d fl h 4u + d u h 4/1 ) + 
+4 > l^u(^—d 4 d 4 h — d a d°h44 — d a d a h — Aa'dji + 3a'd 4 h 44 + 

+d a d T h° T + 2d°d 4 h a4 + Aa'd a h 4 a) ~ 2 V(2(a') 2 + a"<f> + cr>') + 
+3/i447^(4(a') 2 + a") - WpKi + 9„V - 9 4 V) + 29 + 
+^'(d 4 h u + 2d°h a4 - dji + 7a' 'h 44 - 2%f - 8J/cr') + 

+l,u ( W - 2f" - 2d a d°f + 2 fa" + 



+8/(a') 2 + 2f/(0') 2 -2^ 



Qa'f 



2. /i4-component 



4>d 4 (d^h - d u h lxv ) + (j)d u (d u h^ - d^h u4 ) - 3o-'09 A1 /i 4 4+ 



+2d^d 4 f - 2a'dJ - <P'd,h 44 + 2^d,j 



3. 44-component 



/x4 



Aa'dji - I2f(a') 2 + 



1 2cf) , (d fi h fl4 + Aa'h 44 - \d 4 h) - 2d^f - 

0. 



+I2<f)h 44 (a , , - r K „ , ^ 2 
-8a' f - U^'ff + n#f - ~Mfh 44 - £/ 



4. equation for the field / 



d ii {d v h liV - d^hu - d^h + 2d 4 h^ 4 ) - d 4 d 4 h + 
+<9^ ju4 (10<7' - 2|0') + d 4 h 44 (4a' - f 0') + d 4 h(-5a' + 



(15) 



(16) 



(17) 
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iKy) + \%6{y - L)) + 



d 2 \ 2 



O 2 -8>y-2^"-0-^(y)- 



^(y - L)) + /' - 2^0') + 2%f" + 2^9 M / = 0, 



where h = 7mjv^ MJV , /i = r y llv h^ v . Note that these equations can also be obtained by linearizing 
the Einstein equations and the equation for the scalar field, which follow from action ([2]), in 
the background defined by solution of form (j3J), (jlj). 

Let us now discuss the gauge invariance of the linearized theory. It is not difficult to check 
that the quadratic action is invariant under gauge transformations of the form 

/ W = f-fa, 

where Vm denotes the covariant derivative with respect to the background metric jmn, if 
functions £ M satisfy the orbifold symmetry conditions 

e(x, -y) = e(x, y), £ 4 (x, -y) = -£ 4 (x, y) (19) 

(here (/) denotes the transformed field). The existence of these gauge transformations is a conse- 
quence of the invariance of action (j2J) under the general coordinate transformations. Analogous 
gauge transformations were discussed in the case of Randall-Sundrum model without the stabi- 
lizing scalar field and in the case of stabilized Randall-Sundrum model [U El El EE] • One can use 
them to isolate the physical degrees of freedom of the fields Iimn and /. Though this problem 
can be simplified. Indeed, for the case of the minimal coupling of the stabilizing scalar field 
there were found gauge conditions, which allows one to isolate the physical degrees of freedom 
in the general case [6]. Using the conformal transformations, supplemented by the coordinate 
transformations \8\, one can obtain from the gauge conditions, corresponding to the minimal 
coupling of the scalar field, the gauge conditions, corresponding to action ([2D and background 
metric Q: 



(> + U) e 2 > 2 / 3 ] =l(u + l)e^f, (20) 
= 0. 



Analogously we can find the substitution, which allows us to diagonalize equations of motion 
(fl5l) - (fr8j) . as well as the second variation Lagrangian: 

1 / 

h^ u = b^ u - -^^hu - 7ao/t> (21) 

where b^ v is a traceless-transverse field. 

In gauge fl2"U]) and with substitution fl2T|) equation ffTol) fulfills automatically. Equation ffTol) 
takes the form 

\ (W V + |i*v) " *V (2(-0 2 + °" + C^j + ^ = 0. (22) 
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Equation for 44-component (fTTj) simplifies considerably, if one rewrites it in the interval 
(0, L) using a new function g = e 2cj ^(j) 2 ^{h^{x,y) + |4) and taking into account the relation 
between the potential V and functions <r, <fi 



In terms of function g, substitution ( 1211 and gauge conditions (1201) take the form: 

V = Kv-\l^~ 2a <P~ 2,Z 3-\lJ^ ( 24 ) 
^44 = e^r 2/3 9-l^ (25) 

' - K w+ 9^ (26) 

V = 0, 6 = 7^ = 0, 0«V = O. (27) 

Note, that the field g is a superposition of component of the metric fluctuations and the 
fluctuation / of the stabilizing scalar field. 

Substitution of (J2"l"]) - (j2"7|) into equation (TT%|) results in the equation, which can be obtained 
by differentiating ( |23l) by y, and in the boundary conditions on the branes: 

r ~ I (Z7473) ^ + (TO) K - 4-')) <?' " = 0, (28) 

(f + I (^73) ^ + (37T4) K - 4ff ')) 9' - d^g\ y=L . = 0. 

For the case of background solution @, ffTUj) the boundary conditions simplify considerably 
and take the form 

l(^P 2 9' + d^g\ y=+0 = 0, (29) 

1(^73)^' - fi/jl^w = 0. 

It should be noted that such a simplification of the boundary conditions takes place for a class 
of background solutions in five- dimensional Brans-Dicke theory, namely, if the equations of 
motion for the background configuration of the fields can be reduced to first order differential 
equations [8]. 



4 Mass spectrum of Kaluza-Klein modes and 
four-dimensional effective Lagrangian 

First let us consider the tensor modes of the field b^ u , which satisfies equation (I2"2l) . To find the 
mass spectrum and wave functions in the extra dimension, let us represent the field b^ u as 

00 

b^(x,y) = J2^ u (x)My), D W = m 2 n b; u (x), (30) 

n=0 



6 



where □ = rj^ u d^d v . Substituting this into fl22l) . we obtain 

e- 2a m 2 J, n + < - U{p'f + la" + 2a'^- ) Vn + ^ = 0- (31) 



The boundary conditions following from this equation can be obtained by integrating ( 1311) in 
an infinitely small vicinity of the points y = 0, y = L and have the form 

2<rty«l w =+o = 0, (32) 
i/4 - 2aVn|j/=L-o = 0. (33) 



It follows from the general theory [TT] that all eigenvalues of the problem under consideration 
are real and positive. Thus, the tensor sector does not contain tachyons. The eigenfunctions 
corresponding to different eigenvalues are orthogonal with the weight function defined by equa- 
tion fl3~Tl) . The eigenfunctions can be normalized as follows 



L 

2a 



<pe a i> n i> k dy = 5 nk . (34) 

L 



The solution for the zero mode is 



Vi VT - f - 



and if u <C k and kL 35 

^b(L) « ^e" fci (36) 

V 

For the case of massive modes equation (|3~T|) in the interval (0, L) can be solved in a standard 
way (see, for example, [U E]) by passing to variable z = ^e ky ~ kL , which leads to 

il> n {z) = z a (AJ a (z) + BJ_ a (z)) , (37) 

where J a (z) is the Bessel function, a = ^,a = 2 + ^|. In the next section we will show that 
one should take kL ~ 35 to obtain the (weak) Newtonian gravity on the brane at L retaining a 
strong five-dimensional gravity. In this case one can use the approximation z\ y= o = ^-e~ kL ~ 
with a good accuracy, which allows one to drop the singular term J- a (z) in ip n (z), because 
B/A ~ e ~ 2akL anc l the corresponding corrections are negligible. The boundary condition at L 
gives us the mass spectrum of tensor Kaluza-Klein modes, which is defined by 

Ja-l (^) = °" ( 38 ) 

Note that in the limit u = we reproduce the equation for the mass spectrum of tensor modes in 
the unstabilized Randall- Sundrum model J\ (■^ 2 -) = [3]. Normalization constant A is defined 
by formula (|34p . and the normalized wave functions of massive tensor modes look as follows: 

i,M =ivr tm (39) 
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At the point y = L we get a simple formula 

ML) = (40) 

Now let us turn to the scalar sector. To find the mass spectrum of the scalar modes defined 
by equation (I2"2"j) we represent g as 

g(x,y) =e ipx g n (y), p 2 = 

As a result equation ( 1231) and boundary conditions ( 1291) for g n (y) take the form 

rf^f^-^H^M)*"-*.-* (41) 

i-G&mfta'n + e~ 2 ^ 2 n g n \y= + o = 0, (42) 



4 (u>+4/3) 



- e- 2 >^ n | y=i _ = o. (43) 



Equation (1411) is written in the interval (0, L), but it can be combined with the boundary 
conditions (i4"2l) . ( fl3l) . which results in a single equation on the circle S 1 : 



A5/3 e -2a\' 2 e" 2,J 
dr-pT- ) " 9 + 4 ) + ( 44 ) 



e" 4,J /8(3cj + 4) 8(3w + 4) 



e 

/ 

With the help of this equation one can show that 

d f ( ^gfc + g (3w + 4) -^9*91,) = (45) 

for n ^ A; (we suppose that fi n ^ 

/ L dy \C~^ Q n + 9 (3W + 4) J = (46) 



., , L 5 / 3 e" 4CT /8(3cu + 4) r . . 8(3u; + 4)_ , , ., 



From (j46j) it follows that the scalar sector has no zero mode with /i = 0, because if fiQ = 0, 
then go(y) =0. If f3f 2 > 0, then /z 2 > 0, i.e. the scalar sector does not contain tachyons. Using 
the formulas presented above, we will normalize the wave functions as 

2 f L . 5 / 3 e" 4<T /8(3o; + 4) r/ , 8(3o; + 4) . . rN \ 

8(3tu + 4) 



27 

Five-dimensional scalar field g(x,y) can be represented as a series 



-5, 



nk- 



y) = Yl ( fn(x)g n (y), (4$ 



n=l 



where four-dimensional scalar fields (p n (x) have the masses fi n . 

With © and ffTUj) equation (l4"T]) in the interval (0, L) takes the form 



g'n + g'n (2* 



u z UJ + 



4 



g n + e^ kL ^ n g r , 



0. 



(49) 



It can be solved by passing to the variable z = B± e k y- kL , After some calculations (which 
are absolutely equivalent to those carried out in [6J for the scalar sector of stabilized Randall- 
Sundrum model), we obtain 



with 



g n (z) = z q (A n J 7 (z) + B n J^(z)) 

6a; + 7 



(50) 



Q 



6k 



'2m 2 



7 



6oo + 6' 

2u + 3 



2u + 2 



Substituting fl50l) into the boundary condition at zero, we can drop the singular term ~ J_ 7 (z) 
(as for the tensor sector), i.e. B n = 0. The boundary condition at L defines the mass spectrum 
and takes the form 



1 + 



J. 



7-1 



k J 4/i n (3u; + 4)' 



/4(2^ + 2)J 
The wave functions g n (y) take the form 

9n{y) = A n {^e k \y\~ kL ) q J,{^e k \y\- kL ) 



(51) 



(52) 



A, 



3/x n J 7 ^ k 
3 

vi /3 k(3u + 4) 



(53) 



7 2 fc 2 

/4 



1 



4(3u + 4)/x r 
3w 2/ 3 2 




where the normalization coefficient A n is derived from ( 1471) . 

To get an effective four-dimensional Lagrangian of the theory (which also allows one to check 
the correctness of normalization conditions fl34"j) and (|47p by checking the coefficients in front 
of the four- dimensional kinetic terms of the fields), we need the second variation Lagrangian of 
the model. The necessary part of the second variation Lagrangian is 



-7 



_ ( _^ x [ eq . <^} uu _ h« x [eq. ([I7j)] + / X [eq. 



(54) 



where we have taken into account = 0, and [eq. ([15]) ] . [eq. ([17]) ]. [eq. (TT81) ] are the left 
hand sides of equations ([T5]) . ([T7]) . ([TS]) respectively. This result is not surprising: in fact, 
formula ( ]54"I) follows from the definition of the second variation Lagrangian. Thus, in principle, 
one can obtain the second variation Lagrangian using only linearized equations of motion for 
the fields. 
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Substituting (130]) and (|48|) into second variation Lagrangian flUj) , taking into account (|24|) - 
( 1271) . (13~4"|) . fj4"Tj) and integrating over the extra dimension, we get 

S *ff = -4 E / ^ (^'^^ + m^ftj,) - 2 E / ^ X + A*W*) • (55) 

fc=0 ^ k=l 

Thus, we have obtained the effective Lagrangian of the theory, which is the sum of the standard 
four-dimensional Lagrangians for the scalar and tensor fields. We note that the kinetic terms 
in the effective Lagrangian have the proper sign, i.e. there are no phantom fields in the four- 
dimensional effective theory. 



5 Interaction with matter 

Interaction of four- dimensional fields b™(x) and tp n (x) with the Standard Model fields on the 
branes is described by the interaction of the fluctuations of five-dimensional gravitational field 
hfa, with matted on the branes, which has the standard form 

h^ u (x, 0)T^yJ ' -det^ pa (ti)d A x + \ J ^(x, L)T^yJ-det^ pa (L)d 4 x, (56) 

and Tf% being energy-momentum tensors of matter on brane 1 and brane 2 respectively. 
We restrict ourselves to matter on the brane at y = L, which is supposed to be "our" brane. 
Substituting expansions (JHOJ) and (HE} into (EEj), taking into account pE]) . (EZ51). (p6l) and (115]) . 
we get the formula describing interaction of tensor and scalar modes with matter on the brane 
at y — L in Galilean coordinates on that brane 



00 



\ I [ ML)bl u (x)T^ + ML)b;^) T ^~ (57) 

- Ib z \ n=l 



V 

2 

71=1 

The coefficient in front of the zero tensor mode b^ v (x) defines the four-dimensional Planck 
mass on the brane, and with the use of (13"6"1) we obtain 

M Pl = $-\L) « ^e kL . (58) 

If ~ 1 TeV" 1 , whereas kL rs 35, then Mp\ ~ 10 19 GeV. Thus, the hierarchy problem of 

gravitational interaction on the brane at y — L is solved analogously to that in the Randall- 
Sundrum model. The coupling constants to matter on the brane have the form 



j>n(L) / A; j 

— ~y^r 1TeV ' (59) 

where we have used formula (HO]) . When u > 35, one gets a ~ 2 and the mass spectrum of the 
tensor modes is approximately the same as that in the unstabilized Randall-Sundrum model 
0j. In this case mi »s 3.8k and can be of the order of 3 — 4 TeV. 
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Now let us turn to the scalar sector. We would like to mention an interesting feature of the 
coupling constants: in principle it is possible that for appropriate values of the parameters 



UV2P0 

for some j. In this case j-th mode of the scalar field does not interact with matter on the brane. 

Now let us estimate the mass and the coupling constant of the lightest scalar mode, - the 
radion. To simplify the analysis we use the "stiff brane potential" limit - /?| — > 00. Let us 
suppose that fii < k. Expanding the Bessel functions in floTl) into a series, retaining the terms 
up to the second order in ^ and solving the resulting quadratic equation, we get 

2k . . 

Hi « -7=. 60 

\/UJ 



Since uj > 35, then fii < k, which confirms the validity of the expansions of the Bessel functions. 

Note, that if /3| — > 00, then g' n (L) = 0, which follows from the boundary condition ( 1431) . 
Thereby, only the first one of the two terms in ( 1571) . describing the interaction of the scalar 
modes with matter on the brane, remains, and the coupling constants look like 

e« = ~Vz 2/3 g n (L). (61) 

Taking into account (1601) and uj > 35, for the lightest mode - the radion - this constant is 
simply 

1 / k 1 nr , . 

ei "-2vi^r~iviv (62) 

which is ei ~ 1 TeV~ l for k 3 ~ v 2 ~ 1 TeV 3 . The radion mass in this case can be of the order 
of hundreds of Ge V. 



6 Conclusion 

In the present paper we discussed a stabilized brane world model in five-dimensional Brans- 
Dicke theory and found equations of motion for the fields describing excitations above the 
background solution. A convenient gauge and a substitution were found allowing one to diag- 
onalize the equations of motion and to isolate the physical degrees of freedom of the model. 
Analogously to the case of the stabilized model with the minimal coupling of scalar field to 
gravity, the tensor sector decouples from the scalar one. For the background solution ([9]) we 
found the mass spectra of tensor and scalar modes and the coupling constants to matter on 
the brane at y = L, where our world is assumed to be. It was shown that the effective four- 
dimensional Lagrangian does not contain tachyons and phantom fields. It turned out that, 
contrary to the stabilized Randall- Sundrum model, the linearized equations of motion for the 
case of background solution ([9]) can be solved analytically for all physically interesting values 
of the parameters of the model, i.e. this model is exactly solvable in the linear approxima- 
tion. This fact can be useful for obtaining estimates of the influence of the extra dimension on 
processes on the brane for different values of parameters of the five- dimensional theory. 

It was shown that for a certain choice of the model parameters the radion mass may be 
of the order of hundreds of GeV, the inverse size of the extra dimension and the masses of 
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tensor excitations being of the order of TeV. The coupling constants of massive tensor and 
scalar modes appear to be of the order of TeV~ l , whereas the coupling constant of the massless 
graviton appears to be ~ , i.e. the hierarchy problem of gravitational interaction is solved 
on the brane at L. 

Finally we would like to note that the explicit form of the background solution for functions 
a(y) (f)(y) was used only for calculating the mass spectrum and coupling constants. All the 
results related to the gauge choice, the diagonalization of linearized equations of motion and 
the structure of tensor and scalar sectors are valid for any scalar field potential in a stabilized 
brane world model in five-dimensional Brans-Dicke theory with background solution of the form 
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